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Ifosoerical Differentiation With Interpolatiaai (DOI or 

SADDI) 
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C&osed 

given x, successive values ©f (h /rl) f* r; (x) 

for arbitrary r are f®«M, a® vail as tha intar* 
##lated value of f(x). 
105 

(2d 4-5) (r + l) + 2d(l-r) milliseconds* where 
d s depth of interj>olation 
r = number of derivatives desired (see-Kote 
2). 
Depends on the aesh siz»e and the depth of 
interpolation « 

Prior to reading in this routine, memory locations 
3 and k must contain the following ? 
3* 00 F 00 tF 

t * initial location for tessporary 

storage 
00 T 00 eF 
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e - difference "between successive 
interpolates (see Hote l)* 



through 15 pins [{r + 1) (2d - r)/2] + 1 

locatio&s "feeginninig at t* 

When the routine is left, f(x), hf x '(%), 

(fa 2 /2t) f^ 2 *(x) ,..., (h r /rl) f^(x) are in 
mmory locations t, t ♦ lj t + 2, ..., t + r, 

respectively* 

With the argument x in A and the first word of 

this routine in mes&ory location %> entry is »ade 

"by* 



•2- 



» 


50 rF 

50 fP 


£ +1 


26 *J 
00 *F 


f + 2 


00 dF 



-where n » number of sueitttervals ove* the iatervml 
[0,1] (see Description}. 
a «= location of the first tabulated value 
f(x A ). 
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This routine employs Neville's interpolation formula 
and its derivatives, and confutes the interpolate and 
derivatives of f(x) from one or more tables stored 
in the memory. 

The domain of f must he scaled so that Q < x < 1* 
The mesh, therefore, will consistof the n + 1 equally 
spaced points Q ~ x < x. < x p < ... < x « 1, 
and the table must consist of the n + 1 function 
values f (x,) at these points, -where - 1 < f(x.) < 1. 
0iven any argument x, subject only to < x < 1 (thus 
we may have x * x for < i < a), this routine 
selects an e<|ual number of points on each side of 
x "whenever possible . For example, if n « 15 and 
d - 5 (thus requiring 5 points), then for 
(a) xax A the routine uses 



(b) 

(c) 



V *i j 



t , Xy and Xj^ 



8 < x < x Q the routine uses 

X 6* ^7' X 3* ^Q* ^^ ^Hq* 
< x < X-a the routine uses 

*U' *12* *& x Xk * ad *13* 



fc l^ 



eternally, this routine earploys no scaling devices 
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sisce all intermediate calculations Must »e la raage 

prevised a fine eaough mesh is takea. The ©ceurraaea 
of overflow, therefore, indicates the mente should 
"fce refined. 



1, 3.4-1, • . . ,1+3-1,1+ J 

f(x,, x 1+1 * •••> x i+4-i> X i+i^ " t3aen 5r ajvillft,s aethad 

©f successive linear intarpolatioix is given lay 
k*i+j " X} f i,i+l,...,i+j-l ~ (x i " x) f i+l,i+2,,..,i+J 

f - ■ SS rWW% m i.j .i . i .ii nw wiMHMit-r-n-.i l -- ,i ,, . ■ n . -l . r- . --!- . i mn. i n, . ., „n jm.n i n. t i in , r, i> , na i, .. , . , , , m .,. m l « > Ml If I I ■ ■> II , 

1,1+1, «*«,l+j 

i+J i 

tk 
aa* the r derivative of this form is 8 

f (r) s ' x i+j " *■' f i,i+l.,.**,d+^l " ' x i * x ^ f i+l,i+2,...,i+j 

1,1+1,. ♦.,!+.} ~ "' ■ *'- ' i »^ - l . i .n l i r-. 1 , .. iW , -» ,.,. - , , „ „„„ , - „ , .„-,„„,■, .r, i , i „ , -. , ' ■„.. , - , i ,„,r m-,;.,^ , 

x, „. -'x. 
1+j 1 



+ 



r {r^l) _ (r-1) 

14-1,1+2,.* ,,1+4 1,1+1, »..,i+J-l 



x •* x 
x i+j 1 



Kfeere both expressions are applied successively 
with j « 1,2,..., d. Replacing x. . - x, = jh "by j/ia 

l+u 1 

aad x ky x^ + x - x. where x. < x < x. - and 

rearran&iag terBS, these eauatlo&e "beeoaet 



f. = f. 

iji+lj o o o ji+j JL^i+lj o o o ^i+j - 1 

+ (l/j)[(i-k) (f -f X 

jt i>i+l^ • " *,i+j-l~ i+l,i+2,... i i+j)] 
(h r /rl) f. (r ? . . . = f (r > 

H- (l/j) [(I - k ) (fg^...^ - fS,i + 2,...,i +J ) 

^r iji+ly.-ji+j-l " i+l J i+2,.., > i+j' 

-h (f (r " l) - f (r " l) )] 
i+lji+2, . ° . , i+j i , i+1, o . . _,i+j ~l y 

which are the forms used by this routine., 

NOTES (l) The routine will be automatically left when- 

ever two successive interpolates differ by less 

-19 
than e < 2 . If it is desired to go to fixed 

depth, set e = 0. 

(2) r and d must satisfy the inequality 
r < d - lo 

(3) If nothing more than the interpolate is 
desired _, set r = 0. 
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46 ?4x, 
L5 2S, 
L4 31 s 
LO 10F 
42 TIL 
L4 10QL 

00 zm 

46 7IL 

F5 55E. 

40 55L 
22 57L 
15 14? 
L0 XOQL 

40 x4f 
32 55s* 

I£ 22L . 

46 5% 

L5 12F 

L0 10©L 

40 12? 
32 64L 
F5 54L 
42 7<& 
L5 2% 
42 76L 

41 155' 
LI 5F 
40 l4y 

L5 y 

L0 If, 
40 UJP 
L3 UF 
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LOCATION 
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69 

70 
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72 

75 

75 
76 

77 
78 

79 
80 
81 

82 

83 
81* 



iA ieeL 

32 6% 

40 lJF 
50 11F 
71 7F 

26 7U, 
L4 F 
L0 P 
50 14F 
74 11F 
66 F 
S5 B?L 
L4? 

F5 15F 

40 15F 

L0 4F 

L4 101L 
40 66l 
L5 TIL 
L4 101L 
40 7XL 
L5 7^L 
lA IOIL 
40 7^L 
F5 1% 
22 6% 
L5 3F 
LO 10F 
40 3F 
F5 10F 
40 XOF 



!Test difference of successive 
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■•interpolate or differentiate 
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I£ 3J 


lo xom, 


ho IF 
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Are »*tce**slve interpolates 
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exit* 
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